Quantum criticality of the disordered topological insulator 
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We study quantum criticality of the disordered topological insulator in two spatial dimensions, 
using a Wilson-Dirac Hamiltonian on a square lattice combined with Rashba spin-orbit coupling 
(SOC). Employing scaling analysis of the localization length, we estimate the critical exponent of 
metal-insulator transitions inherent to the system to be v ~ 2.7, which is consistent with that 
of the symplectic class. The phase diagram classifying different disordered topological phases is 
deduced from the size-scaling. In the phase diagram, the symplectic metal phase predominates 
over the weakly disordered region due to the effects of carrier-doping and Rashba SOC. In the 
generic carrier-doped case, i.e., away from the particle- hole symmetric point, our data implies a 
quantitative disagreement with a simple interpretation of the disorder-induced topological phase 
in terms of renormalization of the effective mass gap proposed earlier. Our model Hamiltonian 
determined by symmetry considerations is optimal for the description of systems under Ci v point 
group symmetry, e.g., of the one implemented in the HgTe quantum well. 

PACS numbers: 72.15.Rn,73.20.Fz,71.70.Ej,73.43.-f 



I. INTRODUCTION 

The concept of the topological insulator has been 
widely accepted in the condensed-matter community 
since the appearance of the Z2 topological insulator 
(TI) fi~— The TI, which can be regarded as a synonym of 
the quantum spin Hall effect in two spatial dimensions, is 
a time-reversal invariant counterpart of the quantum Hall 
effect. The quantum Hall effect requires a finite, rather a 
strong magnetic field that breaks the time reversal sym- 
metry. The quantum spin Hall effect occurs without a 
magnetic field, induced solely by the spin-orbit coupling 
(SOC) that preserves time-reversal symmetry— The TI 
has been first experimentally observed in a HgTe/CdTe 
quantum well^ after a theoretical prediction of Bernevig, 
Hughes and Zhang (BHZ)il Recently, a InAs quantum 
well has been observed to be a Tl4~— 

Robustness against disorder is a defining property of 
the topological quantum phenomena. Quantifying the 
role of disorder has also played a central role in the con- 
ceptual development of the topological insulator. Ran- 
dom matrices are classified into three categories: or- 
thogonal, unitary and symplectic, according to their 
symmetry^ The distinction between quantum Hall ef- 
fect and quantum spin Hall effect might be most accen- 
tuated in this context. The quantum Hall effect, breaking 
time-reversal symmetry, belongs to the unitary symme- 
try class, while quantum spin Hall systems arc symplectic 
since SOC is activated. A more exhaustive classification 
of the disordered systems a la Refs. [ID and was then 
employed to a complete classification of topological in- 
sulators and superconductors.— ~— QSHs belong to class 
"AH" in this classification. 

It is naturally an interesting question whether criti- 



cal properties of a disordered topological insulator sug- 
gests the existence of a new universality class reflecting 
the nontrivial topological nature of the clean system. So 
far researches in this direction have been performed on 
the Kane-Meleji2 and Z2 network^ - — models. Detailed 
analyses in Refs. I20M22I imply that the critical behavior 
of the disordered topological insulator is undistinguish- 
able from that of the conventional symplectic systems.— 
Furthermore, recently, localization in three-dimensional 
topological systems including weak TI a 24 ' 25 and topolog- 
ical crystalline insulators 2 ^ have been studied. 

As for the phase structure, the disorder-induced TI 
has attracted attentions. The topological non-triviality 
of TI stems from an inverted band gap A due to the 
strong SOC. It has been pointed out recently that disor- 
der can drive this band inversion converting the system 
from an ordinary insulator to a "topological Anderson in- 
sulator (TAI)" .— The nature of this TAI phase has been 
already much discussed^ - — For example, "Is TAI a dis- 
tinct phase?"— has been one of the main questions. We 
have partly addressed this question in Ref. [29| from the 
viewpoint of phase diagram in the (A, IF)-plane, where 
IF is a measure of the strength of disorder. There, related 
to the above issue, we have predicted a series of multiple 
transitions involving TAI. The idea of the TAI has been 
equally applied to the three dimensional systems £2r— 

Here, let us mention the role of Rashba-type SOC, 
which is associated with the structural inversion asym- 
metry of the quantum well. The Rashba SOC is also said 
to be "s z non-conserving" , mixing the decoupled spin j" 
(s z = +1) and I (s z = —1) sectors. We study the critical 
behavior of the disordered TI by focusing on the phase 
boundary between the TI and a metallic phases. How- 
ever, in two-dimensional disordered TI do not exhibit a 
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metallic phase, if Rashba-typc SOC is not existing^ 7 - - — 
In the presence of Rashba SOC, carrier-doping effect also 
plays an important role on the phase diagram since a 
metallic phase is allowed. 

In this paper, we employ a BHZ-like effective Hamil- 
tonian implemented on a square lattice, reinforced with 
the Rashba-type s z non-conserving SOC. Randomness 
is introduced by on-site disorder potentials. We have 
carefully studied the scaling behavior of the localization 
length, using the transfer matrix method, standard and 
believed to be most accurate. A number of different 
disorder-induced quantum phase transitions are identi- 
fied by the analysis of this data. From this analysis, 
we deduce the critical exponent characterizing the dis- 
ordered TI to metal transition of a HgTe quantum well. 
We pursue our analysis to obtain different cross sections 
of the phase diagram, revealing effects of Rashba SOC 
and carrier-doping. Moreover, we revisit the nature of 
TAI in the light of the mass renormalization, focusing on 
the carrier-doping effect. 

The paper is organized as follows. In Sec. II, a model 
Hamiltonian of a HgTe quantum well is shown. The de- 
tailed derivation of this Hamiltonian from the symmetry 
consideration is shown in Appendix In Sec. Ill, crit- 
ical exponent of the metal-insulator transition, which is 
determined by the finite-size scaling, is discussed. Fur- 
thermore, effects of carrier-doping and Rashba SOC on 
the phase diagram is clarified in Sec. IV. Finally, we 
summarize our results in Sec. V. 



II. MODEL 

We start with introducing the following BHZ-like ef- 
fective Hamiltonian supplemented by spin-flip hopping 
terms represented by r(fc), which stems from Rashba- 
type SOC, 

/Mfc) F(fc) \ 
H[k) - (^rt(fc) h*(-k)J- (lj 

The basis is taken to be (|+ f), |— t): 1+ I), \— D), where 
± and f (i) denote the sign of parity and z-component of 
spin, respectively. Spin-up part h(k) of the Hamiltonian 
is given byi 

,,,, fA-Bk 2 -iAk+ \ 

h ( k )={ lAk _ -A + Bk 2 )> W 

with k = [k 2 + k 2 ) 1 / 2 and k± = k x ± ik y . The spin-flip 
part r(fc) is given by 

r <»> = Cf - A) - (3) 

with the coupling constant a of Rashba SOC. 

In the actual calculation, we use a tight-binding model 
corresponding to the above Hamiltonian, substituting 



ki — > sinfci and k 2 — > 2(1 — cosfe^). The tight-binding 
version takes the following form: 

H = 2J [cUrCr + (clt x C r+a + c\.t y C r+b + ll.C.)] , (4) 
r 

where r is a position on a square lattice, and a = (1, 0), 
b = (0, 1) are the primitive translational vector. The 
lattice constant is set to be unity. The derivation of on- 
site potential e r and hoppings t x and t v in the present 
model is shown in Appendix [XJ 

III. CRITIC ALITY 

In this section, we discuss the critical exponent of 
metal-insulator transition in the disordered BHZ model, 
where Rashba SOC is taken into account. We show the 
extensive data on the critical exponents, and compare our 
results with those on other related models in the light of 
universality class arguments. 

In the disordered system, on-site energy e r fluctuates 

as 

_ (A — 4B + W+ \ , . 

€r -y o -b + AB + W-J' lj 

where is a probability variable obeying the uniform 
distribution which takes a value within [-W/2, W/2]. W 
corresponds to the strength of disorder. In the numerical 
computation, the parameters are set as A = B = 1 and 
D = a' = a' n = 0. 



A. Localization length and finite-size scaling 

We briefly review our method to determine a criti- 
cal point and the corresponding exponent from finite- 
size scaling42r— In disordered quasi-one dimensional sys- 
tems, wave functions ipi( x ) ,s decay as ipi(x) ~ e ±zAi ; 
where Ai (Ai < ••■ < Al) is the localization length 
of the system. The largest localization length Al di- 
verges at the critical point of metal-insulator transition 
as Al ~ \q — 9c| _I/ , where q is a parameter representing 
disorder strength W or energy E, and q c is the critical 
value. The largest localization length divided by the sys- 
tem width Al/L, where L is the width of the system, be- 
comes scale- free near the critical point. Therefore, Al/L 
is expanded as 

. N M 

1 = A o + E a «(« ~ 1c) n L n ^ + J2 b n(l ~ q c ) n L n,v+y , 

(6) 

The third term in the above expression is an irrelevant 
length scale collection with a negative irrelevant exponent 
y. In L — > oo, the single-parameter scaling recovers: 

A N 

I ^A / = Ao + 5>n(g-<Zc)"£" A '. (7) 
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Metal-Z 2 


Metal-OI 


A 


1.824 ±0.018 


1.825 ±0.007 


W c 


5.932 ± 0.006 


9.780 ± 0.008 


V 


2.77 ±0.23 


2.60 ± 0.06 


y 


-1.57 ±0.37 


-3.21 ±0.74 


N d 


69 


59 


X? 


0.88 


1.02 



TABLE I. Criticality of metal-topological insulator (Metal- 
Z2) and metal-ordinary insulator (Metal-OI) transitions. 
Critical values of localization length divided by the system 
size Ao, disorder strength W c , the value of critical exponent 
v, irrelevant exponent y, degrees of freedom TV<j, and reduced 
chi-square \r- 



The present system has 2gL localization lengths, where 
g = 4 is the number of internal degrees of freedom at 
each site. These are calculated by the transfer-matrix 
method (See Appendix |C|) . The numerical data of the 
largest localization length is fitted by Eq. ©. As a 
result, the critical value q c and the critical exponent v 
are obtained. In the actual calculation, wc fix the fitting 
parameters to TV = 4 and M = 2. Additionally, some 
remarks on the size-scaling is discussed in Appendix [Dj 



B. Critical exponent 

Now, we turn to discussion of the critical exponent at 
the transition between an ordinary /topological insulator 
and a metal. The largest localization length divided by 
the system width Xl /L and the single-parameter scaling 
part A' in the vicinity of topological (ordinary) insulator- 
metal transition are shown in the upper (lower) panels of 
Fig [1] In this calculation, the standard error of \l/L 
less than 0.015 is adapted. At A = 2 and a — 0.5, the 
system is a topological insulator in the clean limit. As 
one increases disorder strength W, the system goes into 
metallic phase at W ~ 5.93 (the left-upper panel). The 
corresponding single-paramcter-scaling part A' is shown 
in the right-upper panel. From the fitting, the critical 
exponent v of metal-topological insulator transition is 
determined as v ~ 2.77 ± 0.23. As shown in the lower 
panels of Fig. [TJ for the larger disorder strength, the 
system turns into an ordinary insulator at W ~ 9.78. 
Similarly to the previous case, from the fitting, v is de- 
termined as v ~ 2.60±0.06. Both for the two transitions, 
the critical value of localization length Ao and the critical 
exponent v are estimated to be Ao ~ 1.82 and v ~ 2.7, 
which are consistent with those in symplectic systems, 
e.g., SU(2) model^ the Z 2 network models ^"^2^ and 
the Z 2 quantum kicked rotator^ Namely, the topologi- 
cal non-triviality of the present system does not affect the 
quantum criticality of metal-insulator transition. The 
obtained critical values and fitting parameters are sum- 
marized in Table HI 



IV. PHASE DIAGRAM 

In this section, effects of disorder on the phase diagram 
are discussed. We make a close comparison of different 
situations i) with and without carrier-doping and also ii) 
with and without Rashba SOC. 



A. Phase diagram at A = —0.1: Disorder-induced 
topological insulator and metallic phase 

Figure [5] shows the phase diagram at A = —0.1, where 
an ordinary insulator is realized for W = 0. The tran- 
sition line between topological and ordinary insulators 
in the case of a = is denoted by circle symbols (red 
line). For a — 0, a topological insulator is realized in the 
isolated region (2 < W < 9) in (W, 75)-plane, i.e., dis- 
order induces a top olog ical insulator, which is consistent 
with Refs. and |28| Introducing Rashba spin-orbit 
interaction a, a metallic phase can appear due to anti- 
localization i 47 ! 48 Actually, in the strongly disordered re- 
gion (W > 5), metallic phase is allowed around the tran- 
sition line for a = 0. The disorder-induced topological 
insulator region is reduced to 2 < W < 7. In addition 
to this, metallic phase widely appear in the weakly dis- 
ordered region (W < 5), and the transition line connects 
smoothly to the clean limit (E ~ 0.1 and E ~ 8). Note 
that the insulating phase does not connect to the case 
of W = for a = 0, while the metallic phase does for 
a ^ 0. 



B. Phase diagram at E = 0: Identity between a 
topological insulator induced by disorder and that in 
the clean limit 



Phase diagram at E = more precise than that in 
Ref. is shown in Fig. [3J A finite metallic region is 
found, which partitions the two topologically distinct in- 
sulating phases. The inner (week disorder) region is a 
topological insulating phase, and the outer (strong disor- 
der) region is an ordinary insulating phase. In the region 
of -1< A < 0, W < 10 for a = 0, and of -0.5 < A < 0, 
W < 7 for a = 0.5, the disorder-induced topological in- 
sulator is realized. This topological insulator is continu- 
ously connected to that in the clean limit (0 < A < 4 and 
W = 0). Namely, a disorder-induced topological insula- 
tor is essentially the same as that in the clean limit 

The transition line between the ordinary insulator and 
metal is located roughly at W = 10, which is nearly 
equal to the band width and independent of A. Since 
Anderson localization occurs when the energy scale of 
disorder is larger than that of band width, metallic phase 
cannot exist for W > 10 even if the Rashba SOC becomes 
stronger. It is expected that as one increases a from 
a = 0, the metallic region gradually spreads around the 
transition line for a — and it converges roughly to 
W < 10. 
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FIG. 1. Localization length divided by the system size (Xl/L) (left) and the single-parameter-scaling part A' (right) for 
A = 2.0. The upper (lower) panels show the localization length in the vicinity of the transition point from a topological 
(ordinary) insulator to a metal. The horizontal axis corresponds to (W — Wc)! 1 '" oc (L/^) 1 ^ . The dashed lines in the left 
panels denote the critical point determined from the the finite size scaling shown in the right panels. 



The inner transition line located between topological 
insulating and metallic phases, on the other hand, depend 
on not only band width but also detailed band structure, 
i.e., a and A. As one approaches A = 4, where the band 
gap closes at W — 0, the critical disorder strength W c 
decreases and the transition line connects to (A, W) = 
(4,0). Topological insulator phases in < A < 4 and 
4 < A < 8 are separated each other, even in the presence 
of disorder. Z2 topological insulators for < A < 4 and 
4 < A < 8 have the same strong topological number 
v = 1. However, for 4 < A < 8, bound states appear 
at a dislocation ) 4 ^ 50 while it does not appear for < 
A < 4. In this sense, these two topological phases are 
distinguished from each other. 



C. Phase diagram at E — 0.5: Carrier-doping effects 

Next, we discuss carrier-doping effects in the disor- 
dered system. Figure U] shows the phase diagram for 
E = 0.5, which corresponds to the case of Fig. [3] with 
carrier doping. In the absence of SOC (a = 0), Ander- 



son localization always occurs, therefore carrier-doping 
slightly affects to the transition line, as shown in Fig. 0] 
(circle symbols with red line). Although for a = the 
region of disorder-induced topological insulating phase is 
suppressed but survives for A > —0.5, that for a = 0.5 
vanishes. A topological insulator exists in the smaller 
phase region (1.2 < A < 3.2) than that for a = 0. 

Compared with the case of absent carrier density (Fig. 

, phase structure in the weakly disordered region dras- 
tically changes in the presence of SOC. For —0.5 < A < 
1.2, the system is metallic at W = both for a = and 
a = 0.5. In contrast to the case for a = 0, the critical 
line for a = 0.5 continuously connects to the critical point 
(A = -0.5 and A = 1.2) at W = 0. For a = 0, a metal- 
lic phase cannot appear due to localization. Therefore, 
only the sign of A is relevant to the topological number, 
i.e., the transition line connects to A = and W = 0. 
On the other hand, for a 7^ 0, since a metallic phase 
can appear due to anti-localization, the transition line 
connects to the band edge for W = 0. As a result, the 
wider metallic region is realized in the weakly disordered 
regime (W < 6) in the presence of SOC. 
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FIG. 2. (Color online) Phase diagram in (W, _E)-plane at 
A = —0.1. The lower panel shows the enlargement of the up- 
per panel. The error bar denotes ±2e, where e is the standard 
error defined in Appendix [Cl 
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FIG. 3. (Color online) Phase diagram in (A, VF)-plane at 
E = 0. 



D. Renormalization of the mass 

The notion of the "disorder-induced" topological insu- 
lating phase is based on the observation that in Dirac sys- 
tems disorder not only renormalizes the magnitude of the 
mass (i.e., the band gap) but sometimes it even changes 
its signi2£ Namely, in the presence of the intra-orbital 
hopping, whose amplitude is B in the Hamiltonian, dis- 
order can renormalize mass A from negative (ordinary 
insulator) to positive (topological insulator) values. A 
simple and useful method to quantify renormalization of 
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FIG. 4. (Color online) Phase diagram in (A, VF)-plane at 
E = 0.5. 

the mass term is the self-consistent Born approximation 
(SCBA) . 28 ' 51 In this section, however, we remark that 
SCBA for the carrier-doped case does not quantitatively 
become a good approximation. 

The self-energy in the disordered system is given by 

within SCBA, where G k (E) = (E-H (k)-'E(E)+i0)- 1 
is the Green's function, and Ho(k) is the Hamiltonian in 
the clean system. The self-energy E(-E) is decomposed 
into T,(E) = T, (E) + T, z (E)a z due to the symmetry re- 
quirement (See Appendix [E]) . The Renormalized mass A 
and Fermi level E are given by 

A(E) = A + RcZ z (E), (9) 
E(E) = E-ReE (E). (10) 

In the renormalization picture, a topological insulator is 
realized when A(E) is positive and E(E) is located in 
the band gap. 

In the following, we discuss the phase boundary in the 
carrier-doped (E = 0.5) and undoped (E = 0) cases. 
The band inversion point A(E) = and the band edge 
E(E) = Eg, where 2E g is the magnitude of the gap, are 
denoted by the (red and green, respectively) solid lines in 
Fig. [5l In the case of E = 0, (the left panel for a = and 
the center panel for a = 0.5), where there is no carrier 
in the system, E is fixed to E = due to the chiral sym- 
metry, therefore the critical point between the ordinary 
and the disorder-induced topological insulators is deter- 
mined only by A = 0. In the less carrier-doped case, 
SCBA gives quantitatively good result in weakly disor- 
dered region (W < 5), even in the presence of spin-orbit 
coupling, as shown in the left and center panels of Fig. [5] 
On the other hand, in the carrier-doped case of E = 0.5, 
SCBA cannot quantitatively give the critical point of the 
topological phase transition. The lines of A = (the 
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FIG. 5. Comparison between the critical points calculated by the size-scaling and by the SCBA. The band inversion point 
(A = 0) and the band edge (E — E s ) calculated by the SCBA are denoted by the (red and green, respectively) lines. The 
critical points determined by the size-scaling are denoted by the circle symbols, which are the same as in Fig. [3] and [4] 



red solid line in the right panel) does not coincide with 
the actual critical points (the circle symbols in the right 
panel) determined by the size-scaling. This mismatch 
suggests that Anderson localization, which is beyond the 
range of SCBA, plays an important role for the topo- 
logical quantum phase transition for carrier-doped cases. 
For the carrier-undoped case, the system is already insu- 
lating in the clean limit and does not directly related to 
the localization. In such cases, SCBA gives the correct 
critical points. 



V. SUMMARY 

In this paper, we have clarified the criticality of a dis- 
ordered two-dimensional topological insulator. The ob- 
tained critical exponent is v ~ 2.7, which is consistent 
with that in symplectic class (All). Namely, topologi- 
cal non-triviality docs not affect the criticality of phase 
transition. The phase diagram of the disordered sys- 
tem has also been studied in the light of the disorder- 
induced phase and effect of Rashba SOC. Carrier-doping 
and Rashba SOC yield a wide metallic region in the 
phase diagram, due to the anti-localization. In addi- 
tion, carrier-doping affects the disorder-induced topolog- 
ical phase, i.e., the renormalized mass A(i?) does not 
quantitatively give the critical points between ordinary 
and disorder-induced topological insulators. 

A two-dimensional topological insulator implemented 
in a quantum well has advantages for experiments. The 
Fermi level and strength of Rashba SOC can be tuned 
by the gating. Moreover, it is possible to tune Rashba 
SOC with tuning structural-inversion-asymmetry of the 
quantum well. We believe that our results obtained in 
the present paper will be verified in such quantum wells 
with HgTe£ and InAs^ii 
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Appendix A: Rashba SOC in a two-dimensional 
topological insulator 

Hamiltonian of HgTc quantum well with inversion sym- 
metry has been derived by Bernevig, Hughes, and Zhang 
(BHZ model) X In this Appendix, we derive Hamiltonian 
with Rashba SOC, from symmetry consideration in C^ v 
symmetry, i.e., on a square lattice with a electric field 
perpendicular to the xy-plane. The tight-binding Hamil- 
tonian of BHZ model^ is given by 

H = ^ [ c UrCr + {cUxCr+a + c\.t y C r+b + h.C.)] , 
r 

(Al) 

where r is a position on a square lattice, and a, b are 
the primitive translational vectors. The internal degrees 
of freedom of the system is labeled by z-component of 
total angular momentum j z = ±1/2 for s-orbital and 
j z = ±3/2 for X ± iF-orbital. The basis is taken as 
(|l/2), |3/2), | - 1/2), | - 3/2)). The system has time- 
reversal O, four- fold rotational R4, and mirror M y sym- 
metries. The corresponding (anti-)unitary matrices act- 
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ing on the internal degrees of freedom arc given by 

9 = -iSyK, (A2) 



a = 



Ra 



My — Sy, 



(A3) 
(A4) 



where a and Sj are Pauli matrices in the orbital and 
spin spaces, respectively. The on-site potential e r should 



satisfy the following relations. 




QcrQ^ 1 = e r , 


(A5) 




(A6) 


M y e r Ml = e r . 


(A7) 


This restricts the form of e r to 




e r = (A - 4B)cT g s , 


(A8) 



where the parameters A and B are real numbers, Simi- 
larly, t x has the following relations, 







(A9) 


R^t x R^ 


= ft 


(A10) 


Myt x Ml 


"X 1 


(All) 



with R\ = <j z s z . Thus, the form of t x is determined as 



A a 
t x = Da s + Ba z s - i—a-ySo + i—<J z s y 



JyCToSy + — (JyS v , 



(A12) 



where all the parameters A, D, and a are real numbers. 
a, a', and qo denote strengths of Rashba SOCs. t y is 
obtained by the rotation of t x as 



ty — R^txR 



L)a so + Bcr z s - i—a x s z - i—a s x 



a 

i-^v z s x 



ao 



(A13) 



Here, we confirm that Rashba SOC docs not appear in 
inversion-symmetric systems. If we consider D^h symme- 



try, two-fold rotation i?2 



s x along the x-axis 



also becomes a symmetric operation. Consequently, the 
following relation holds. 



R2t x Rl — t x 



(A14) 



This leads to a = a 1 = ao = 0. To be sure, Rashba SOC 
is not allowed in inversion-symmetric systems. 

In the momentum space, the Hamiltonian with Rashba 
SOC reads 

H(k) = (A - 4B)a z s 

+ 2(£>ct + Ba z )so(cos k x + cos k y ) 

+ [Aa y s - (aa z + a'a )s y ] smk x 

+ [Aa x s z + (aa + a'a z )s x ] sink y 

-)- OLQOySy (cos k x — cos k y ) . ( Al 5) 
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FIG. 6. (Color online) Energy spectra of the lower conduction 
band with a = (left) and a = 0.5 (right). The parameters 
are given as D = a' — ao = 0, A = B = 1, and a = 0.5. 



In the above expression, the linear terms of k near T- 
point with spin flip \—{aa z + a'ao)s y smk x + (aero + 
ct'u z )sx sin k y ] are present, differently from the case in 
Ref. l53l This is because the symmetry in a square lattice 
(Ci V ) is lower than that in k p theory (Coot,). Actually, 
in the axial symmetric case (Coav), the spin- flip term for 
j z = ±3/2 states satisfies 

(k, 3/2|JT|fc, -3/2) = e m (k', 3/2\H\k', -3/2), (A16) 

where k 1 = e~ % i* k. As a result, the leading term is 

(fe, S/2\H\k, -3/2) oc ki, (A17) 

namely, the linear term is prohibited. The linear term is 
possible in the discrete symmetric case C^ v . 

Appendix B: Clean limit 

In this Appendix, we clarify the electronic properties 
of a topological insulator with Rashba SOC in the clean 
limit. Hamiltonian of the system is given by Eq. (|A15[) . 

1. Electronic states in the clean limit 

Figure [5] shows energy spectra of the conduction band 
for a — and a — 0.5 in W = 0. For a = 0, the bottom 
of the conduction band is located at T point. When a 
is switched on, spin-degeneracy is lifted and the minima 
of the energy band become located at finite fc-points as 
shown in the right panel of Fig. [5] 

The corresponding density of states is shown in Fig. 
[71 which is shown only for E > since it is symmetric 
with respect to E — 0, due to the chiral symmetry. The 
behavior is consistent with that in Ref [48|, where the s- 
orbital is focused on and the p-orbitals are integrated out. 
The conduction band exists about in 0.1 < E < 8. In 
the presence of Rashba SOC, the density of states takes 
a larger value at the band edge {E ~ 0.1 and E ~ 8), 
due to the multiple-minima of the conduction band. In 
the continuum limit, the density of states at the band 
edge diverges since the shape of energy band becomes 
wine-bottle. 
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FIG. 7. (Color online) Density of states in the pres- 
ence/absence of Rashba SOC a. 



2. Phase diagram in the clean limit 

Before discussing the disorder effects, we clarify the 
phase diagram in the clean limit. When s z is con- 
served (a = 0), spin Hall conductance is quantized. It 
is known that such a topologically non-trivial phase is 
robust against weak Rashba SOC, which breaks spin ro- 
tational symmetry, i.e., s z is not conserved. In this case, 
although the spin Hall conductance is not quantized, we 
can define the Z 2 topological invariant as discussed by 
Kane, Mele and Fui^ However, as one increases Rashba 
SOC, the magnitude of the band gap decreases and van- 
ishes at the critical value a = a c . For a > a c , the system 
becomes metallic and irrelevant for topological quantum 
phenomena. Figure [5] shows the phase diagram deter- 
mined by evaluating the magnitude of the band gap. In 
the absence of a,— ^ a Z 2 topological insulator is realized 
in < A < AB with a xy = e/27r and AB < A < 8B with 
<7 xy = — e/27r, where a xy is a quantized spin Hall conduc- 
tance with respect to z-component of spin. At A = 0, AB, 
and 8B, the band gap vanishes (zero-gap semiconductor) 
and the phase transition occurs. An ordinary insulator 
is realized in A < and A > 8B. In the presence of a, 
the gapless points A = 0,4B, and 8B still remain, and 
a Z 2 topological phase can exist about in E < 0.5 for 
A = B = 1 and a = 0.5, as shown in Fig. 13 where the 
transition line between Metallic and insulating phases is 
denoted by the solid line. If one further increases a, the 
region of Z 2 topological phase shrinks, and vanishes at 
a > a c = A. Moreover, the phase diagram is symmetric 
with respect to A = 4 due to the chiral symmetry in the 
case with D = a' = qq = (See Appendix IE]) . 



FIG. 8. (Color online) Energy (E) - mass (A) phase diagram 
in the clean system. Z2 and OI denote Z2 topological insulator 
and ordinary insulator respectively. 



of our numerical simulation. In addition, we show the 
definition of the standard error. 

Let us consider the following Schrodingcr equation dis- 
cretized on a two-dimensional lattice: 

Eipi,j = ei,jipij + t x ipi+i,j + 4^/-l,J 

+ t v i>j,j+i + tl^i,j- 1 . (ci) 

This can be rearranged into a one-dimensional form in 
terms of the transfer matrix Mj as 




(C2) 



where ipj = (ipi.i, ■ • ■ iPi,l) T - Mj is generally a 2gL- 
dimensional matrix with g being the number of internal 
degrees of freedom (g = 4 in the present model), and is 
given by 



Mr 



t x tyl+ t x tyl — 







(C3) 



where ej = diag(e/,i, • • • ,ej,i), {1+)i,j = I- = 

it . Due to the disorder, wave function decays exponen- 
tially as ipi ~ e ±/ / Ai , where is a localization length. 
Xi is deduced from the transfer matrix as 



A, 
L 



In 



,(* = !,••• ,2gL) 



(C4) 



Appendix C: Transfer- matrix and the numerical 
method 



The transfer matrix method is used for estimating nu- 
merically the localization length of the system. Here, we 
describe some details of the method, and the protocol 



where to,; is the eigenvalue of Af( Lx ' = Ml x _i ■ ■ ■ Mi. 

Thus, to estimate A, we have only to calculate the 
eigenvalues of the product of transfer matrices, M^ Lx ' = 
ri/=i Mr- A direct calculation, however, usually fails, 
since numerical error tends to accumulate in the product 
of matrices. A common resort to this failure is the use 
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of Q-R-dccomposition, i.e., after an iterative use of the 
decomposition as 

M x = Q x Rx, (C5) 
M I+1 Q I = Q I+1 R I+1 , (C6) 

one finds 

= Af ix • • • M x = Ql„R Lx ■■■Ri, (C7) 

where Q/'s and i?/'s are respectively unitary and trigo- 
nal matrices. Neglecting Ql x in the above equation turns 
out to be a fairly good approximation for large enough 
L x . Then, A/( L =) can be regarded product of trig- 
onal matrices. This means that one can safely estimate 
the eigenvalues of this matrix simply by multiplying the 
corresponding diagonal element of each trigonal matrix 
Rj. Recall that the eigenvalues of a trigonal matrix are 
its diagonal elements, and a product of such trigonal ma- 
trices are also trigonal. 

The actual calculation has been done by following the 
steps as listed belows^ 

i. Set the initial value of the unitary matrix Q, e.g., to 
be 1, for simplicity. Set also as d = (0, • • • , 0), e = 
(0,--- ,0),L X = 0. 

ii. Generate the rii transfer matrices and calculate the 
product M = M ni ■ ■ ■ M X Q. 

hi. QR-decomposition. M = QR. 

iv. Store the diagonal components of R as d a — d a + 
In \r a \,e a = e Q + (ln |r a |) 2 , where r a is the a-th diag- 
onal component. Calculate the system length L x : 
L x = L x + n i . 

v. The localization length and the approximated stan- 
dard error e a are given by X a /L = L x /(d a L) and 
e a = L(X a /L) 2 y/e a /Ll - (d a / L x ) 2 n t / L x respec- 
tively. 

vi. If the error e a is larger than the accuracy goal eo, 
go to step ii. 

In step iv. we have introduced a parameter rii . This must 
be chosen appropriately, in accordance with the nature 
of the problem. If rii is too large, the resulting numerical 
error would become intolerable. An appropriate choice 
of rii can be made by tuning d a to be symmetric, i.e., d a 
should appear as a pair such that d a = — d b . This comes 
from the fact that the eigenvalues also appear as pairs 
e ±L x /Ai ; . . . ^ e ±ix/A L _ Here, we have chosen as rii = 5. 
and have checked that the symmetry of d a = —d b is sat- 
isfied as (d a + d b )/(d a - d b ) < 0.02. 

Last but not the least, we have also employed a parallel 
computation. Using A^cpu = 64 CPUs, we calculate the 
localization length of -/Vcpu independent systems with 
the accuracy of ~ e x (Ncjpu) 1 / 2 in each calculat- 
ing node (i = 1, ■ ■ ■ ,Ncpu)- Collecting the results as 



t _ v^JVcpu T (i) j _ r^fcpu _ x^Nq-pv Ai) 

we obtain the average Xl/L and the error . This paral- 
lclization is actually not that trivial, since the j-th "sam- 
pling" does depend on (j — l)-th sampling by the unitary 
matrix Q. Yet, we believe that the parallelization is fully 
justified for sufficiently large L l x /rii. 



Appendix D: Remarks on the choice of boundary 
conditions: determining the critical point 

In the actual calculation, finite-size effect can become 
serious problem to determine the critical point. For an 
insulator-insulator transition or for the weakly disordered 
case, it is useful to apply an open boundary condition, 
in order to reduce the finite-size effect, as we shall see in 
the following. 



Periodic boundary condition 

Two panels of Fig. [3] represent a typical example of 
our size-scaling analysis. The left panel shows the local- 
ization length (divided by L) as a function of disorder 
strength as varying the size L of the system. To perform 
numerical estimate of the localization length we make 
the system in the form of a long, quasi-one dimensional 
tube of length L x ~ oo and circumference L (periodic 
boundary condition is applied in the circumferential direc- 
tion). Here, the Rashba SOC is finite (a = 0.5), and the 
phase diagram shows a metallic region between the two 
- one topological, the other ordinary - insulating phases 
(the system belongs to the symplectic symmetry class). 
Furthermore, in the range of W relevant in these plots, 
the system undergoes a quantum phase transition from a 
metal to an ordinary insulator. This appears, in the left 
panel of Fig. |Hl as the change of the size-dependence of 
Xl/L from an increasing to a decreasing behavior, with 
a size-independent critical point in the middle, located 
roughly around W ~ 9.8. More precise determination of 
the critical point needs, however, single-parameter scal- 
ing analysis as shown in the right panel of Fig. [5J due 
to a non- negligible finite-size effect. Dashed lines in the 
left panel indicate the position of such a critical point: 
(W c , A ) = (9.82 ± 0.20, 1.76 ± 0.17), determined by the 
size-scaling. 



Open boundary condition for insulator-insulator transition 

Next we consider the case of absent Rashba SOC: 
a = 0, i.e., the system belongs to the unitary class. In 
this case, only direct transitions between the two distinct 
topological insulators are expected; no metallic phase ap- 
pears in between. This is problematic to our numer- 
ical analyses, since fitting the data to identify a criti- 
cal point is much more difficult between two insulating 
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phases. Imagine that on both sides of the scale invari- 
ant critical point Xl/L is a decreasing function of the 
size as far as the system bears no surface states. To 
overcome this difficulty, we adopted, here in FIG. [TQ1 an 
open boundary condition, since then the system can sup- 
port a pair of gapless edge states in the quantum spin 
Hall phase, i.e., it becomes a metal- insulator transition. 
Fitting data to determine a critical point is much eas- 
ier in the vicinity of a metal-insulator transition than in 
insulator-insulator transitions ! 5 ' The left panel of FIG. 
1101 demonstrates determination of such a transition from 
a one-dimensional metal to an ordinary insulator (E = 0, 
A = 0). 

One might notice, however, presence of a stronger 
finite-size effect here compared to the case of a ^ (FIG. 
GO left panel). Here, the existence of a scale invariant 
point is no longer evident, and besides, apparent loca- 
tion of the focal point is much deviated from that of the 
critical point: (W C ,A ) = (7.461 ± 0.046, 2.137 ± 0.053) 
(indicated by dashed lines) determined by the single- 
parameter scaling (fitting shown in the right panel). Ex- 
istence of such a stronger finite-size effect is a disadvan- 
tage of an open boundary condition! 45 ' 56 



Open boundary condition for weakly disordered case 

Use of an open boundary condition has another area 
of utility. In the weak disorder region, the localization 
length Xl becomes large. The larger Xl becomes, the 
larger system one needs to consider to perform estimates 
of the localization length. Use of an open boundary con- 
dition is useful for such a case, since the typical length 
scale of the system, i.e., the critical value of the localiza- 
tion length Ao is smaller in system with an open bound- 
ary. This is, on the other hand, applicable only for a 
metal to an ordinary insulator transition; e.g., in FIG. [3] 
only the phase boundary between OI and metal could 
be determined in this way. This is because a metal- 
topological insulator transition in an open boundary con- 
dition becomes metal-metal transition, which is also diffi- 
cult to determine the critical point as insulator-insulator 
transition. 



Appendix E: Chiral symmetry 

Symmetry argument is useful to clarify structures of 
energy bands and phase diagrams, without detailed cal- 
culation. In this Appendix, we show that the form of the 
Green's function is restricted by the chiral and particle- 
hole symmetries. Due to these symmetries, energy spec- 
trum should be particle-hole symmetric and phase dia- 
gram should be symmetric with respect to A = AB. 

Wh en D = a 1 = and W = in Eqs. (jJT]), (|A8]l. 
(IA12[) . and (|A13|) . the system has the chiral symmetry 



as 

f HY t = -H, (El) 

rc r r^ = x s y c r . (E2) 

As a result, the energy spectrum becomes particle-hole 
symmetric. 

In addition to this, one can prove that the phase dia- 
gram is symmetric with respect to A = AB for D = a' = 
0. Let us define the local gauge transformation as 

S CjjJ St = {-l) I+J a x s y c IiJ: (E3) 

which corresponds to (7r, 7r)-shift in the momentum space, 
where r = (/, J) is a position on the square lattice. Re- 
placing A — AB — > — (A — 4_B), and operating the above 
transformation, the transformed Hamiltonian is equiva- 
lent to the original one: 

Si?| A _ 4B ^_( A _ 4S )S T = H. (E4) 

Note that although the disorder potential term does not 
microscopically invariant for the above transformation, 
disorder-averaged quantities such as the Green's func- 
tion become invariant. Consequently, the phase dia- 
gram becomes symmetric with respect to A = AB. Note 
that a chiral-asymmetry may affect the behavior of gap- 
less/gapful edge states^ 

Appendix F: Symmetry of the local Green's function 

From the symmetries of the system, the form of the 
local Green's function is restricted as follows. The lo- 
cal Green's function is defined by Gj^j> r (z) = (r, j z \(z — 
H)- 1 ^,^), with j z ,j' z = ±1/2, ±3/2. In the follow- 
ing, this is regard as 4 x 4 matrix G{z). Similarly to 
Eq. (|A8|) , from the four- fold rotational and time- reversal 
symmetries, G{z) is decomposed into 

G(z) = G (z)a a s Q + G z (z)<r z s . (Fl) 

In addition to this, if the system has the chiral sym- 
metry shown in the previous section, the particle-hole 
symmetry is also present. The particle- hole symmetry is 
expressed as 

CcX- 1 = ct x c*. (F2) 

This leads to 

-G(-z*) = a x G*(z)o- x . (F3) 

For z ~ iO, with the help of the above relations, the local 
Green's function satisfy 

RcGo(iO) = 0, lmG z {i0) = 0. (F4) 

Namely, the renormalized Fermi level is given by 

E(0) = Tr Re G(i0) /A = 0. (F5) 
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FIG. 9. The localization length and its single-parameter scaling in the presence of Rashba spin-orbit interaction (a = 0.5) are 
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